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ABSTRACT 

To obtain an unbiased sample of bright Lya blobs (^Lya ^ IC'^^ergss"'^), we have undertaken a 
blind, wide-field, narrow-band imaging survey in the NOAO Deep Wide Field Survey Bootes field 
with the Steward Bok-2.3m telescope. After searching over 4.82 deg^ at z = 2.3, we discover four 
Lya blobs with LLya = 1-6 - 5.3 xlO''^ ergss"^, isophotal areas of 28 - 57n", and broad Lya line 
profiles {Av — 900 - 1250 km s^^). In contrast with the extended Lya halos associated with high-z 
radio galaxies, none of our four blobs are radio-loud. The X-ray luminosities and optical spectra of 
these blobs are diverse. Two blobs (3 and 4) are X-ray-detected with Lx{2—7 keV) = 2-4 xlO^"* 
ergss^^ and have broad optical emission lines (C IV) characteristic of AGN, implying that 50% of 
our sample blobs are associated with strong AGN. The other 50% of blobs (1 and 2) are not X-ray or 
optically-detected as AGN down to similar limits. The number density of the four blobs is '--^ 3 x 10^^ 
Mpc~'^, comparable to that of galaxy clusters at similar redshifts and 3x lower than that found 
in the SSA22 proto-cluster at z = 3.1, even after accounting for the over-density of that region. The 
two X-ray undetected blobs are separated by only 70" (550 kpc) and have almost identical redshifts 
(< 360 kpc along the line-of-sight), suggesting that they are part of the same system. Given the rarity 
of the blobs and our discovery of a close pair, we speculate that blobs occupy the highest density 
regions and thus may be precursors of today's rich cluster galaxies. 

Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium 



1. INTRODUCTION 

The recent discovery of extended Lya nebulae, the so- 
called Lya blobs, may ultimately provide clues to how 
galaxies formed in early Universe (z '--^ 2 - 5). The 
most luminous Lya blobs have been discovered mostly 
by deep narrow-band imaging surveys and have typi- 
cal sizes of ~ 10" (~ 1 00 kpc) and Lya luminosities 
of ~ 10'^^ ergss-^ (e.g iKeel et all [19991 iSteidel et al " 
Francis et al.l[200ll iPalunas et al.ll2004l : iDev et al 
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Smith k JarvisI 120071) . However, what powers the 



blobs and what they will evolve into remains a mystery. 
To begin to answer these questions, we need first to con- 
strain their basic properties, including how common they 
are, how they are spatially distributed, and how diverse 
their possible power sources may be. 

Despite intense interest, these characteristics, for even 
the brightest Lya blobs, are poorly constrained due to 
the ab sence of an approp riate sample. Gian t, luminous 
blobs (jSteidel et al.ll2000l ) appear to be rare (jSaito et al.l 
I2OO6I ) at z = 3 — 5, so an efficient, large- volume sur- 
vey is required to detect any. Furthermore, the survey 
must be blind to be representative of blob statistics and 
must include ancillary data (in the X-ray, other optical 
bands, etc.) to better constrain possible sources of blob 
emission. To date most Lya blobs have been discovered 
by targeting t wo p reviously known ove r -dense regions 
(ISteidel et al.l 12000 : Franci s et all 120011: IPalunas et all 
12004 iMatsuda et al..,2004.) . iMatsuda et all (|2004D dis- 
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covered 33 additional Lya blobs in the proto-cluster re- 
gion at z = 3.1 where the two giant Steidel et al. blobs 
reside, and spectroscopic follow-up confirms that the 
blobs are located in three filamen tary structures trace d 
by compact Lyg-emi tting galaxies (jMatsuda et al.|[2005[) . 
IPalunas et al.l (|2004l ) report that four Lya blobs are asso- 
ciated with an over-dense region of Lya emitters around 
the galaxy cluster J2142-4420 at z = 2.38. The first and 
only study of the environment of a Lya blob f ound with- 
out a ny prior knowledge of its surroundings ()Dev et al.l 
120051 ) shows it to lie in a 3x over-dense regio n traced by 
Lya-emitting galaxies (jPrescott et al.|[2008[) . Although 
these results suggest that bright Lya blobs occupy higher 
density regions, we still lack an unbiased, statistical mea- 
surement of how frequent and clustered they truly are. 

The origin of the blobs' Lya emission is simi- 
larly murky. Proposed mechanisms include grav- 
itational cooling radiation from ac creting gas 
(iHaiman. Spaans. fc Ouataertj l2000l : iFardal et all 
I2OOII), galactic superwind s driven by a starburst 
(Tan iguchi fc Shioval |2000D. photo-ionization by ex- 
tended star f ormation (^Matsuda et al] |2007|), and a 
hidden AGN (|Haiman fc Rccs 2001). ISmith fc Jandi 
()2007t l report the discovery of a giant Lya blob that 
appears to have cooling origin, while others suggest 
that blobs are exp eriencing outflows dWilman et al.l 
l2005h or are static (jVerhamme et al] l2006t ). It is also 
not clear what central sources the blobs harbor, e.g., 
what fraction of those central sources contain AGN? 
These ambiguities arise from the opt ically-thick nature 
of the Lya line (cf. lYang et ahl 120061 ). which clouds the 
interpretation of the blob kinematics, and from the lack 
of multi-wavelength data, which, were it available for a 
representative blob sample, could reveal the nature of 
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any central sources. 

To obtain a unbiased sample of bright Lya blobs ideal 
for determining their properties and, ultimately, for re- 
solving the origin of their extended emission, we are con- 
ducting a blind, wide-field, narrow-band imaging survey 
for the brightest (^Lya ^ 10'*'* ergss"*) and most ex- 
tended {Aiso ^ 50 O") object s , sim ilar to those orig- 
inally found by ISteidel et all (|2000[ ). To date, only 
five s uch blobs have been discovered among several sur- 
veys fSteidel et all '2000: 'FranciseLaD |200H iDev et all 
[2005t ISmith fc Jarvis 2007). We select the NO AO Deep 
Wide Field Survey (N DWFS) Bootes and Cetus fields 
(jJannuzi &: Devlfl999D . a s well as the Cosmi c Evolution 
Survey (COSMOS) field (iScoville et al.l[200l as our tar- 
gets, given the complementary, multi-wavelength data 
available for these fields. We choose z = 2.3 as our sur- 
vey redshift so that the optically thin Ha A6563 emission 
line from the surrounding gas falls in a relatively clean re- 
gion of the NIR sky spectrum. The comparison between 
Ha and Lya, as well as analysis of other lines including 
[O II] A3727, [O III] A5007, and Civ, will enable us in 
future papers to stu dy the kinematics of t he surrounding 
gas in more detail (jPiikstra et al.l [20061 ) and to better 
constrain the nature of any central sources. 

In this paper, we present the initial results from our 
Lya blob survey in the NO AO Deep Wide Field Survey 
(NDWFS) Bootes field. In Sj^l we describe our narrow- 
band imaging survey design, observations, and selection 
of Lya blob candidates. In 331 we present the first re- 
sults from the survey and its spectroscopic follow-up. We 
also discuss the number density and environment of the 
Lya blobs in comparison with previous studies. In 21 we 
summarize the results. Throughout this paper, we adopt 
cosmological parameters: Ho — 70/i7okms~* Mpc~*, 
riivi = 0.3, and f^A = 0.7. All magnitudes are in the 
AB system (IOkelll974n . 

2. OBSERVATIONS 
2.1. Survey Design 

Using the 90P rime one-degree field imager 
([Williams et al .' 2004) mounted on the Steward Observa- 
tory Bok 2.3m telescope on Kitt Peak, we carry out our 
narrow-band imaging survey with a custom narrow-band 
filter (hereafter AfB403 or NB). This narrow-band filter 
has a central wavelength of Ac ~ 4030A, designed 
for selecting Lya-emitting sources at z « 2.3, and a 
band-width of AAfwhm ~ 45A that provides a line 
of sight depth of Az ~ 0.037 corresponding to 46.8 
h^Q Mpc at z = 2.3 in the comoving frame. The 
90Prime camera consists of a mosaic of four 4k x 4k 
blue sensitive CCDs with a plate scale of 0'.'45 per pixel. 
The detector field-of-view is ^ 0.8 deg^ after accounting 
for cosmetic defects in two of the CCDs. This wide 
field-of-view is ideal for our survey as the large area 
coverage compensates for the narrowness of the filter, 
producing a survey comparable in volume to the largest 
intermediate band filter survey (ISaito et al.l[2006h . 

To date, we have finished narrow-band imaging obser- 
vations with a total sky coverage of ~ 12 deg^ in three 
fields where extensive ancillary data sets are publicly 
available: the NOAO De ep Wide Field Survey Bootes 
and Cetus fields (NDWFS: [Jannuzi fc De "v[ ll999ti and the 
Cosmic Evolution Survey field (COSMOS: iScoville et"an 
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Fig. 1. — Sky coverage of our narrow-band imaging survey for 
Lya blobs within the NOAO Deep Wide Field Survey (NDWFS) 
Bootes field. The shaded region represents the area covered in our 
Bok 2.3m -|- JVB403 filter survey, and the dot-dashed lines repre- 
sent the sub-fields of the NDWFS. The gaps in the sky coverage are 
due to the large gaps between the 90Prime CCDs and a large trap 
in one of the CCDs. The total sky coverage is 4.82 deg^, exclud- 
ing bright stars and artifacts, making our survey one of the widest 
FOV, narrow-band imaging surveys to date. The four circles rep- 
resent the Lya blobs that we have discovered (see i|3.2l for details). 
None are radio-loud. The easternmost and westernmost blobs are 
detected in X-rays, while the the pair of blobs in the upper middle 
region is not. 



[200l . In this paper, we present the initial results from 
the survey within the NDWFS Bootes region centered 
at a = 141^ 32'" 05!7, 6 = +34° 16' 47'.'5 (J2000). We 
carried out these observations from April to June 2007 
over 16 dark nights. Figure [1] shows our sky coverage 
overlayed with the NDWFS field layout. After exclud- 
ing regions with artifacts and bright stars, we achieve 
a total sky coverage of 4.82 deg^. Although the expo- 
sure times range from 6 to 20 hours, all the Lya blobs 
in our final sample would have been detected even in the 
shallowest sub-fields ( §2.2p . suggesting that trading solid 
angle for depth will improve constraints on the demo- 
graphics of Lya blobs in future surveys. The sky cov- 
erage corresponds to a comoving survey volume of 2.1 
X lO^hj^Mpc^, which is ^ 16x larger than the area 
cover ed by the deepest nar row-band Lya blob survey to 
date ([Matsuda et al.ll2004[ ). 

We reduce the data with the IR AF mscred mosaic 
data reduction package ([Valdeslll998f ) . The data are cor- 
rected for cross-talk between amplifiers, bias-subtracted, 
and dark-subtracted. For flat-fielding, we use twilight 
flats together with night-sky flats, which are median- 
combined from object frames without the bright stars 
each night. Satellite trails, CCD edges, bad pixels, and 
saturated pixels are masked. The astrometry is cali- 
brated with the USNO-Bl.O catalog ([Monet et al.ll2003^ 
using the IRAF ccmap task. The images are scaled us- 
ing common stars in each frame and stacked to remove 
cosmic rays. Finally, the images are transformed into 
the same world coordinate systems as the NDWFS Bw 
band images, which are resampled with a coarse 90Prime 
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Fig. 2.— The (B^"* - NB403) color-magnitude diagram for all 
sources detected in either the Byf and NB images within the Bootes 
survey field. Right and top axes show the corresponding equivalent 
widths in the observed frame and NB fluxes, respectively. The 
sharp boundary at JVB403 > 24.0 is due to sources detected only 
in the Bw band. We select line-emission objects with the criteria 
NB < 23.5 (dot-dashed line) and B^"' - NB > 1.5 (dashed line; 
^l^obs > 140A). Squares represent the final Lya blob candidates 
(see FigO. 

pixel scale (0'.'45 pixel"^). We observed 3-5 spectro- 
photometric standard stars in each night to derive ex- 
tinction coefficients and zero points for the JVB403 mag- 
nitudes. The major drawback is the poor seeing, which, 
in the final combined images, ranges from I'.'S to 2" with 
a median of I'.'T. However, as we discuss in the see- 
ing does not impact our ability to identify the blobs with 
sizes greater than 250". 

To identify blob candidates requires that we subtract 
the continuum emission within the NBA03 bandpass. We 
estimate the continuum using archival, deep, broad-band 
Bw images from t he NO AO Deep Wide Field Survey 
data set (BwAUK: iJannuzi fc Devlll999() . The wide Bw 
filter (Ac « 4135A and AA « 1278A) encompasses our 
JVB403 filter near the central wavelength. The 5-a lim- 
iting magnitude of the Bw images is ~ 26.6 mag for a 
2" diameter aperture. The seeing of the NDWFS Bw 
images ranges from O'.'S to 1'.'47 with a median of 1'.'2. A 
total of 20 NDWFS tiles cover our narrow-band imaging 
fields as shown Figure [T] 

2.2. Selection of Lya Blob Candidates 

To find Lya blob candidates, we construct photometric 
catalogs in the JVB 403 and Bw bands using SExtractor 
(jBertin fc Arnouti 1996). Because our narrow-band im- 
ages are much shallower and typically have worse seeing 
than the broad-band images, we make "detection" im- 
ages (NB+B^fif) by adding the ArB403 and Bw images 
after scaling them according to their signal-to-noise ra- 
tios (S/N). After identifying sources in the "detection" 
images that have least 5 pixels that are 1.5cr above the lo- 
cal sky, we run SExtractor in double-image mode on the 
AfB403 and Bw images with these detection images. In 



other words, we first find the sources in the detection im- 
ages and then obtain photometry at their position in the 
NB403 and Bw images to make two separate (narrow- 
and broad-band) catalogs. We adopt Kron-like elliptical 
aperture magnitudes (i.e., MAG_AUTO in SExtractor) to 
derive photometric properties. Our use of the "detec- 
tion" images ensures that 1) all the sources detected in 
either the NB or Bw band are included in our catalog 
and 2) the elliptical apertures determined by the more 
extended sources in (JVB+Bw) are large enough to in- 
clude all the light from both the NB and Bw images. 
This last point is critical because we do not attempt to 
match the seeing between the two bands at this stage. 

The selection of Lya blob candidates from the NB and 
Bw photometry catalogs consists of two steps: 1) selec- 
tion for line (hopefully, Lya) emitting objects with large 
line equivalent widths and 2) selection for spatially ex- 
tended objects with a larger angular extent in line emis- 
sion than in the broad-band. 

First, we choose candidates by requiring that they are 
detected above the completeness limits of the JVB403 im- 
ages (JVB403 < 23.5 mag).'^ We also require that candi- 
dates have observed-frame equivalent widths larger than 
140A [EWrcst > 42A) corresponding to B™"* - NB > 
1.5 (vertical and horizontal lines in Fig. [2]). Candi- 
dates must be extended (SExtractor CLASS.STAR > 
0.95 in both the NB and Bw bands). There are a total 
of ~ 450 objects over 4.82 deg^ satisfying these criteria. 
Because most of the Lya blobs (~30/35) discovered by 
iMatsuda et all (|2004D have EWohs > 200A (SWrcst > 
50A), our selection criteria would include them if blob 
properties remain the same between z = 2.3 and 3.1. 

We estimate the continuum flux density, fcontJ ^^'^ 
continuum-subtracted line flux, Fune, of these objects 
using the following relations: 
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(1) 



where Fb and Fnb are the total flux in each filter derived 
from the Bw and JVB403 magnitudes, respectively. A As 
and AA nb represent the band- widths of the Bw and NB 
filters, respectively. Figure [5] shows the (B^"* — NB) 
color excess as a function of JVB magnitude for all ob- 
jects detected in either the JVB and Bw bands within our 
survey area. Here, B^"* represents the broad-band con- 
tinuum magnitude corrected for the emission line: B^"* 
= -2.5 log -48.6.^ 

At our survey redshift, the only possible interlopers 
are nearby [O II] A3727 emitters at z « 0.08. However, 
such objects rarely have equivalent width s larger than 
lOOA in the rest frame (|Hogg et al.lll998[) . Therefore, 
the contamination of our z = 2.3 Lya source catalog by 
nearby star forming galaxies is expected to be minimal. 
Even if we reduce our selection limit to EWohs ^ 65A, 
we end up with the same final Lya blob sample and none 
of our conclusions in this paper change. 

Second, we identify those line-emission selected objects 

^ The completeness limit varies by ±0.5 mag depending on the 
sub-field due to seeing variations and different exposure times. 
However, our final blob sample does not depend on the choice of 
completeness limit because all the blobs are at least ~ 1.5 mag 
brighter than this cut (JSj. 
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Fig. 3. — Distribution of isophotal areas and Lya luminosities of 
Lycf-emitting sources (circles) selected with the first criteria (Fig. 
[2]l. The horizontal and vertical dashed lines represent the selection 
criteria for the final Lya blob candidate sample; Ai^o > 25 □" and 
L(Lya) > 1.5 X 10*^ ergss"^, respectively. Gray dots represent 
the Aiao-LLyo relation for simulated point sources. Note that our 
four final blob candidates (filled circles) are located well above this 
relation (the solid line represents the upper 2cr limit), indicating 
that they are extended sources despite the poor seeing. Simulated 
point sources mix with extended Lya-emitting sources below our 
selection limits on Aiso and Lj^ya, making it diflicult to distinguish 
blobs from point sources there. The squares represent the predicted 
sizes and Lya luminosities of the 35 Matsuda ct al. (2004) blobs 
were they observed with our observational set-up at z = 2.3. Our 
final blob candidate selection criteria (dashe d lines) would find 6 - 
8 (filled squares) of the brightest and largest iMatsuda et al.l II2004I ) 
Lya blobs (see i|3.3l l. 



that are more spatially extended in Lya than their con- 
tinuum counterparts (Fig. [3]).'* We measure the spatial 
extents of the Lya emission in the continuum-subtracted 
images. After registering the JVB403 and Bw images at 
the sub-pixel level and matching their seeing, we con- 
struct continuum-subtracted JVB403 images by applying 
the relations in Eq. ^ in 2-D. We measure the isophotal 
area of the emission region by running SExtractor with a 
threshold of 5 x 10~^* ergss"^ cm^^ arcsec^^. This mea- 
surement threshold is ^ 2.3 x higher than that adopted 
by Matsu da et al. (2004), but because our survey red- 
shift {z = 2.3) is lower than theirs {z — 3.1), we gain 
a factor of ~ 2.4 in surface brightness. Although this 
measurement threshold is comparable to the rms noise 
of the continuum-subtracted NB403 images (~ 0.5 - 1 
xlO~^^ ergss^^ cm~^ arcsec"^ depending on the field), 
we are still able to measure the sizes over this measure- 
ment limit because the bright central parts of the objects 
were already detected with high S/Ns. 

One of the potential problems in detecting a blob is 
the contamination of our extended Lya candidate sam- 
ple by point sources, especially given the poor seeing of 
the narrow-band images. To quantify this effect, we first 
place artificial point-sources with a range of luminosi- 



* ISaito et al.l II2006I ) selected spatially extended objects by re- 
quiring the FWHM in their intermediate-band to be larger than 
that in the broad-band image, but the large seeing difference be- 
tween our JVB403 images (~1'.'7) and the NDWFS Bw images (~ 
1'.'2) prevents us from adopting this approach. 



NDWFSJ14aBp3531 



A 



NDWFSJ1431p3531 



iir 



10 20 30 40 50 60 10 20 30 40 50 60 

Size [arcsec^] 

Fig. 4. — An example of recovery test for the faintest blob 
(blob 2) in our final sample. We show the distribution of the re- 
covered isophotal area for two NDWFS sub-fields for illustrative 
purposes. This blob is actually located in the NDWFSJ1431-I-3531 
sub-field (right) and has the worst recovery fraction (/rccv) in the 
NDWFSJ1428-I-3531 sub-field (left). The bins with zero size indi- 
cate non-detections. The vertical dashed line represents our selec- 
tion criteria for blob size (25n"). 

ties (-^Lya — - 10*^ ergss^^) into the sky regions 
over the whole NDWFS field and measure their sizes and 
fluxes in the same manner as for the extended sources. 
We then determine isophotal area and line luminosity 
limits above which extended and point sources can be 
differentiated. 

Figure [3] shows the distribution of the angular sizes 
and line luminosities of the 185 Lya-emitting candidates 
assuming that they are all located at z = 2.3. The 
open and filled circles represent the line-emitting objects 
selected using the line-emission criteria (Fig. [2]), and 
the gray dots show the relation between the sizes and 
brightnesses for the artificial point sources. As final blob 
candidates, we select objects with isophotal areas larger 
than 25 □" and line luminosities brighter than 1.5 x 10"*^ 
ergss^^. Below these limits, extended and point sources 
mix, and the sizes of blobs cannot be measured reliably 
as explained below. 

Because the chosen isophotal threshold is comparable 
to the rms sky noise, wc test how reliably we can measure 
the spatial extent of the blob candidates. We cut out a 
small (51x51 pixels) section around each blob candidate 
from the {NB — -Bw) image, filter it with a smoothing 
kernel,^ place each postage stamp into 300 - 1000 empty 
sky regions in each of the 20 NDWFS sub-fields in Figure 
[U and extract the sources with SExtractor in the same 
way as for the real data. Then we check how often and 
accurately their sizes are recovered from these simulated 
images. This procedure tests how the detectability of a 
given blob candidate changes across the fields due to the 
seeing and exposure time variations. 

Figure H] shows an example of this recovery test for the 
faintest blob in our final sample (blob 2, see We 
show the distribution of the recovered isophotal sizes of 
the artificial blob for two sub-fields. While most (>95%) 
of time, the artificial blob is detected, its size is mea- 
sured with a large spread because of non-uniform sky 
background noise. We adopt this spread of the recovered 
size distributions as the error in the blob size, Aiso- 

A recovery fraction of a blob (/rocv) in a sub-field is 
defined as the fraction of time that the artificial blob is 
recovered with a size larger than 25 0". These recovery 
fractions are averaged over the whole field to obtain the 



We use a 3 X 3 pixel convolution mask with a FWHM : 
which is a default kernel in SExtractor. 
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Fig. 5. — Images of the four Lya blob candidates. From left to right arc the AfB403, NDWFS Bwi aiid continuum-subtracted (JVB— Bw) 
images. The ticks are spaced every 10", which corresponds to ^ 82 physical kpc at ^ = 2.3. The elliptical apertures used for photometry 
and blob identification are shown in the left panels. The Bw images are smoothed to match the seeing of the JVB images. In the right 
panels, the contours represent the surface brightness limit of 5 X 10~^* ergss"^ cm~^ arcsec"^, within which we determine the blob sizes. 
The location of the 1.5"-wide slit used for spectroscopic follow-up with the 6.5m MMT is shown for blobs 1, 2, and 3 in the Bw image. We 
spectroscopically confirm that all four blobs lie at z = 2.3 (Fig. \E^. Blobs 3 and 4 are broad-line QSOs that are both detected in X-rays 
[Lx(2-7keV) = 1.6 - 4.1 x 10***ergss-i]. None of the blobs are radio-loud. 



final recovery fraction of the blob for this survey. The size 
error becomes comparable to the measured size below 
our sample selection criteria and the recovery fraction 
drops from >90% to ^50% at the selection boundary 
for the brightest blob candidates. The recovery fractions 
are used to correct the incompleteness of our survey in 
calculating the blob number density in 



3. RESULTS AND DISCUSSION 

The selection criteria, EWohs > 140A, Aiso > 25 0", 
and L{Lya) > 1.5x lO'^^ ergs s~^, yield four extended 
Lya candidates at z = 2.3 with sizes of 28 - 57 O", Lya 
luminosities of 1.6 - 5.3 x 10^'^ ergss^^, and EWohs of ~ 
170 - 420A (Table [U. Figure E] shows the JVB403, Bw, 
and continuum-subtracted (JVB— Bw) images for these 



candidates. Except for Blob 4, which appears to have 
more than two continuum counterparts, we are able to 
identify single host galaxies in the Bw bands for the 
other three blob candidates. Although these continuum 
sources are barely resolved in Bw, the line-emission re- 
gions are extended over ~ 5-10" in the (JVB — Bw) im- 
ages. It is noteworthy that all four blob candidates have 
obvious continuum counterparts in the deep broadband 
images. To date, all Lya blobs at z = 2 — 6, except the 
blob discovered bv lNilsson et al.l (2006), have continuum 
counterparts within the extended emission or nearby. We 
also show the locations of the blob candidates on the sky 
in Figure [H Because blobs 1 and 2 are separated by only 
70" (~ 550 physical kpc a.t z = 2.3), they are shown in 
the same panel in Figure [S] 



6 



Yang et al. 



Lya NV 



1.5 

1.0 

0.5 

0.0 
1.0 

0.8 
0.6 
0.4 
0.2 

0.0 
6 
5 
4 
3 
2 
1 

4 




I ' ' ' ' I ' ' ' ' I 
.blob 2 





[0 II] 



H|3 [d :ftl] 



1200 1210 1220 1230 

Rest Wavelength [A] 



4000 



4500 



5000 



5500 



6000 



6500 



Wavelength [A] 



Fig. 6. — Extracted ID spectra for four blobs from 6.5m MMT long-slit spectroscopy (blobs 1, 2, and 4) and the AGN and Galaxy 
Evolution Survey (blob 3; Kochanek et al. in prep.). For blobs 1 and 2, there is no other emission line visible up to 6500A. The broad 
Lya profiles and the absence of other emission lines indicate that blob 1 and 2 are true Lya blobs. These two blobs have almost the 
same redshifts (z = 2 .318) within Sz ~ 0.001. Their spectra do not show any AGN signatures (e.g., ClvA1549). Neither of these blobs is 
X-ray-detected ( i|3.2l l. Blobs 3 and 4 show other broad emission lines (e.g., C IV A1549 and/or C III] A1909), indicating that these blobs 
host QSOs. Unlike blobs 1 and 2, these blobs are X-ray detected [Lx(2— 7keV) = 1.6 — 4.1 X 10** ergss"^]. Because of the blending in 
the Bw images, the blob 4 spectrum contains an interloper, an [O ll] emitter at 2 ~ 0.266 (H7, H/3, and [O in] lines). The shaded regions 
indicate where strong sky lines are present. 



3.1. Spectroscopic Confirmation 

We spectroscopically confirm tlie redshifts of all four 
blob candidates using spectra taken with the 6.5m "JVIul- 
tiple Mirror" Telescope (JMJVIT) Blue Channel Spectro- 
graph and an existing optical spectrum from the AGN 
and Gal axy Evolution Su rvey (AGES; Kochanek et al. 
m prep., ICool et aW^m h. For blobs 1, 2 and 4, we ob- 
tained low-resolution long-slit spectra on UT 2008 Febru- 
ary 13, April 29, and June 7. We used the 500 gpm grat- 
ing with a 1'.'5 X 180" slit, which provides a dispersion 
of I.2A per pixel and a spectral resolution of 5.4A over 
the wavelength range 3400 ~ 6400A (1000 - 1900A in the 
rest frame) . Because blobs 1 and 2 are separated by only 
~ 70", we observed them simultaneously for a total of 
6 hours with the slit position (P.A. = —12.1°) shown in 
Figure \5[ Blob 4 is the most ambiguous object in our 
sample as its continuum counterpart is not clear due to 
the blending in the image. Therefore, we put a slit 
at P.A. ~ 90.0° to include all the continuum components 
seen in the Bw band images. The spectra were reduced, 
extracted, and flux-calibrated in the standard manner 
with I RAF. 

The MJVIT spectra confirm that blobs 1, 2 and 4 are 
indeed z = 2.3 Lya blobs. We show the extracted spec- 
tra of these three blobs in Figured In each case, there is 
an emission line at ~ 4030A, which agrees well with the 
central wavelength of the JVB403 filter. For blobs 1 and 
2, there are no other emission lines over the entire wave- 
length range, which covers the redshifted [O II] A3727, 



H/3 A4868, and [O III] A5007 lines from any z = 0.08 
star-forming galaxies. The Lya lines have large velocity 
widths of ~ 900 and 1250 km s~^, and observed-frame 
equivalent widths of ~ 150A and 190A, respectively. The 
broad line profile and the absence of other emission fea- 
tures in blobs 1 and 2 indicate that the emission line is 
Lya. The stellar continuum of blobs 1 and 2 bears a 
strong resemb lance to the UV sp ectra of Lyman break 
galaxies (e.g.. IShaplev et aT1l2003f l. Notably, there is no 
evidence of broad emission at the wavelengths of Si IV 
and C IV, so the presence of an unobscured AGN is ruled 
out with confidence. A detailed analysis of the contin- 
uum and spectral features of these blobs will be presented 
in a future paper. For blob 4, we identify two continuum 
sources, an [O II] emitter at z ~ 0.266 and a QSO that 
produces broad Lya and C IV lines at z ~ 2.32. The 
Lya line has a velocity width of ^ 950 km s""'^ and an 
observed-frame equivalent width of ^ 85 A. Note that the 
EWs measured from the long-slit spectra are lower than 
those estimated from the narrow-band imaging due to 
slit-loss. 

For blob 3, we obtain a spectrum from AGES, which 
has extensive optical spectroscopy for all the X-ray de- 
tected galaxies with Jab ^ 22.0 mag in the Bootes field. 
This spectrum shows that blob 3 is at z = 2.32 and 
contains a QSO with strong, broad Lya emission (Fig- 
ure [S]). Blob 3 is perhaps similar to the diffuse Lya ha- 
los associated with some QSOs Ce.g.. lBunker"erall[2Q0l 
iWeidinger eiralll20f)5h . 
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3.2. Multi- Wavelength Properties of Confirmed Lya 

Blobs 

Are any of the four blobs associated with strong ra- 
dio or X-ray s ources? None is d etected in the VLA 
FIRST survey (iWhite et al.l [l99l at 1.4 GHz at the 
detection limit of 1 mJy. Assuming a power-law spec- 
tral energy distribution of S{h') cx v^^'^, this limit cor- 
responds to a rest-frame 1.4 GHz luminosity density of 
3.2x10^^ ergss^^Hz^^, which is ~ two orders of magni- 
tude fainter than the powerful radio galax ies associated 
with Lya halos (e.g.. iReuland et al]l2003f l. However, in 
X-rays, the properties of these blobs vary wildly. Blobs 
3 and 4, whose optical spectra have broad QSO lines, are 
detected in the Chandr a XBootes survey (jKenter et al.l 
[20051: iBrand et al.l[200fih with L(2-7keV) = 1.6 - 4.1 x 
10^'* ergs s~^. Blobs 1 and 2 are not detected at a similar 
depth and have upper limits of L(2-7keV) < 0.63 and 
0.18 X 10^** ergs s~^, respectively. Therefore, at least two 
of the blobs (50%) are associated with AGN. Table[T]Hsts 
the optical and X-ray properties of the four blobs. 

It is unknown what powers the copious Lya emission 
of blobs. Previous observational studies suggest diverse 
mechanisms, e.g., a cent ral A GN (Dev et al.l2005l) . galac- 
tic superwi nd (e.g.. iWilman et al. 20051. extended star 
formation (iMatsuda et al.l 2007f). and cool i ng radiation 
(iNilsson et al.l l2006t ISmith fc Jarvii l2007t ISmith et all 
|2008() . vet there is no smoking gun. Due to its blind 
strategy, large volume, follow-up optical spectroscopy, 
and overlap with wide Chandra X-ray imaging coverage, 
our work here is the first to place unbiased limits on the 
fraction of blobs with luminous AGN (50%). The di- 
versity of the blob optical and X-ray properties is also 
interesting, suggesting that AGN like those detected in 
half the blobs are not necessary to power the others. 

3.3. Rarity of Lya Blobs 

Our survey indicates that extended Lya nebulae are 
extremely rare at a redshift of 2.3. Based on the com- 
pleteness (recovery) test in H2.21 we estimate the num- 
ber of Lya blobs within the survey volume as N = 
Si l//recvj where f}^^^ is the recovery fraction (Table[T]). 
The discovery of only four Lya blobs over 4.82 deg^ yields 
a number density of 2. 5 ±1.1 x 10~^ /ifo Mpc~^ for blobs 
with Aiso > 25 □" and L(Lya) > 1.5 x lO''^ ergss"!. 
Note that this number density is comparable to those 
of galaxy clusters in the nearby a nd high-z Un i verse , 
n ~ 10-5 - 1 0-^ h^Q Mpc-^l'e.g.. iBahcall eTall [20031 : 
lPapovichll2008D . Because our survey fails to find blobs as 
bright (LLya > 10^4 ergss-1) and l arge (A;,„ > 150 □") 
as the brightest and largest in the ISteidel et al.l (j2000l ) 
sample (thei r Blob s 1 and 2), we conclude that the 
ISteidel et al.l (|200iol) blobs are even rarer objects with a 
number density of n < 0.5 x 10^^ h^f^Mpc'"^. 

To compare ou r blob number density with those from 
previous surveys ( Matsuda et al.l i2004: Saito ct al. 2006), 
we determine how many of the IMatsuda et al.l (|2004l ) 
blobs would be detected in our survey if they were lo- 
cated at z = 2.3. Using the co ntinuum subtra c ted N B^^^^ 
images of the 35 blobs from IMatsuda et all (|2004l ). we 
scale the blob surface brightnesses and sizes to z = 2.3 
according to the adopted cosmology and assuming that 
their physical sizes and Lya luminosities do not change 
from z = 3.1 to 2.3. These images are convolved with 
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Fig. 7. — Number density of Lya blobs at different redsliifts. 
Open and filled squares represent the number density from our 
narrow-band imaging survey at z = 2.3 for the blobs with and 
without X-ray detections, respecti vely. We also show tw o number 
density estimates for the blobs in IMatsuda et ahl I I2004I) . with or 
without correction for the over-density of the SSA22 proto-cluster 
region. We incorporate the same luminosity-size limits on Matsuda 
et al.'s result as obtained in this paper for a proper comparison. 
For blobs at hig her redshifts, we show only u pper (3 < z < 5 ; 
ISaito et al] [20061) and lower limits (z = 6.595: fOuchi et aLllMOSl) 
from the Lyo blob searches in the Subaru/XMM-Newton Deep 
Survey field (SXDF). 

Gaussian kernels to match our poor seeing, rebinned to 
the QOPrime pixel scale, and given Poisson noise. We do 
not account for the difference between the filter (JVB403 
vs. NBA97) band-widths, because the resolving power of 
two filters is similar. We place the simulated images into 
our continuum-subtracted iVB403 images and measure 
their sizes and luminosities in the same way as described 
in ^12.21 We repeat this process 500-1000 times to derive 
the range of recovered luminosities and sizes. 

We show the distribution of sizes and Lya luminosi- 
ties {Aiso - Li^ya) in Figure [31 Based on this test, we 
expect that 6 - 8 of the 35 IMatsuda et all (|2004D blobs'^ 
are detectable given our luminosity-size criteria, yield- 
ing an effective blob number density of 4.6-6.2 xlO^^ 
/lypMpc^'^. Considering the recovery rate of each Mat- 
suda blob, which is -100% for i(Lya) > 5x lO'*^ 
drops to 50% at our selection boundary, we expect to 
detect on average ~ 6.0 of their blobs with our survey 
criteria. Once we factor in the difference between our 
survey volume and Matsuda et al.'s, we should have de- 
tected —98 blobs — not four — provided that the blob 
number density at z = 3.1 and 2.3, and between the two 
survey fields, is the same. 

The comparison of our results with those of lSaito et al.l 
(|200a is more difficult because we do not have high 
S/N measurements of their blob surface brightness pro- 
files. Therefore, we estimate the blob number density 
for L(Lya) > 1.5 x lO^^ergss"^ from the brightest bin 
in their Figure 13. We consider the resulting number 
density of ~ 6.7 x 10~^ Hjq Mpc^'' an upper limit at 
z = 3-5, given that some fraction of their blobs might 
not have been detected with our survey criteria. 

We show the number density of Lya blobs from the dif- 
ferent samples and at different redshifts in Figure [3 after 
applying the same luminosity-size limits as our survey. 

^ In the detectabilit y simulation, the avera ge values of extracted 
Aiso and LLya for six IMatsuda et al] II2004I ) blobs (their LAB 1, 
2, 3, 4, 8, 10) satisfy our selection criteria. Two other blobs (their 
LAB 6 and 15) are located on our selection boundary, but satisfy 
the selection criteria more than 50% of the time in the simulation. 
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Open and filled squares represent the number density es- 
timates from our survey with and without the two X-ray 
detected blobs. For the Matsuda et al. sample, we show 
both the measured number density and that scaled down 
to compens ate for the known over-density of the SSA22 
field {5 ~ 6: ISteidel et"ani2000( l. For Saito et al.'s survey, 
we show the upper-limit discussed above and the red- 
shift range (horizontal error bar). We also show a lower 
limit on the number density at z — 6.595 derived from 
the recent discovery of a L ya blob in the Subaru/XMM- 
Newton Deep Survey field ijOuchi et al.|[2008 ). 

At face value, Lyg blobs within th e z = 3.1 proto- 
cluster from the iMatsuda et al.l ()2004f ) survey are ~ 20 
- 30 times more numerous than blobs in our blind field 
survey. Even if we account for the galaxy over-density in 
the proto-cluster, the blob number density in our survey 
is still by a factor of 3 lower, suggesting that extended 
Lya emission is closely related to (or enhanced by) a 
clustered environment. However, we cannot rule out that 
this discrepancy might reflect real evolution between z = 
2.3 and 3.1. 

3.4. A Close Pair of Lya Blobs 

One of our most interesting results is the discovery 
of the pair of blobs (1 and 2) separated by only 70", 
which corresponds to ^550 physical kpc at z — 2.3. 
To estimate how unlikely it is to find two blobs within 
70" of each other, we calculate the expected number of 
galaxy pairs assuming that the spatial distribution of 
Lya blobs can be represented by a two-point correla- 
tion function, £,{r, z) — [r/ro(z)]~''', where rQ{z) is the 
scale length of galaxy clustering at the redshif t z. By 
projecting this correlation function onto the sky (|LimbeH 
Il953l ). we derive the angular correlation function uj{9), 
which measures the excess probability above random of 
finding a galaxy at an angle 6 from another galaxy. For 
a given power-la w correlation function, one can derive 
uj{9) = ^^6li-T (jPeeblesI [l98fll : lEfstathiou et al.l [l99lh 
such that 



A,,, — 



rl{z)d'a-''{z) 



dN 



dz 



H{z) 



dz 



dN , ■ 
——dz 

dz 



(2) 

where dc[z) is comoving distance at z, dN/dz is the num- 
ber of galaxies per unit redshift interval, H{z) is the Hub- 
ble parameter at z, and = 0rr[(7 — l)/2]/r(7/2). 
We adopt dN/dz as a top-hat function for our survey 
redshift interval and assume a fiducial value of 7 — 1.8. 

From this angular correlation function, we estimate the 
number of galaxy pairs within Op using: 

■'m , , ,,2Tiede 

-N{N-l)[l+u:{e)]——, (3) 



Np{< dp) 



where fl is the survey area, and N is the number of 
blobs found in our survey. For the typical galaxy cor- 
relation length, rp 7 Mpc in the comoving frame (e.g., 
iMaddox et al.lll990D . we predict a negligible pair count, 
Np ~ 0.005. The expected number of close pairs does 
not increase dramatically as the correlation length in- 
creases. Even for the correlation l ength of the richest 
galaxy clusters , tq ~ 30 Mpc (e.g.. fBahcall et aP 120031 : 
lPapovichl l2008) for rcdshifts out to z ~ 1.5, only Np ~ 
0.05 pairs are expected. Therefore, we conclude that the 
observed close pair of Lya blobs is unlikely to occur for 



a reasonable range of clustering strengths if we assume 
that the two Lya blobs belong to individual dark matter 
halos. 

Increasing the size of our blob sample (Yang et al. 
2008, in prep.) will better constrain the clustering of 
Lya blobs. For the time being, it is intriguing that 
we do discover a close pair. Our MMT spectra con- 
firm that the two blobs have almost the same redshifts 
within Sz ~ 0.001 and a corresponding line-of-sight sep- 
aration of Sr ~ 360 physical kpc assuming no peculiar 
velocity. The similarity of their redshifts, as well as 
their small separation on the sky (corresponding to 550 
physical kpc), leads us to speculate that these two blobs 
might lie within a single dark matter halo, i.e., within 
a massive proto-group or cluster of galaxies. Note that 
the analytic ACDM model predicts a virial radius (r2oo) 
of 430 kpc at z = 2.3 for a dark m atter halo with a 
mass of 10" Mq (IMo fc White! [200I) . Previous studies 
also suggest that L ya blobs are phenomena that favor 
dense environments (Matsuda et al." 2005t iPrescott et al.l 
2008) . This hypothesis could be tested by using a deeper 
narrow-band imaging survey to characterize the environ- 
ment of this blob pair with the spatial distribution of 
fainter, smaller blobs and/or faint, compact Lya emit- 
ters. 

4. CONCLUSION 

In this paper, we present initial results from our 
blind, wide-field, narrow-band imaging survey in the 
NOAO Deep Wide Field Survey Bootes field to constrain 
the number density, environment, and multi- wavelength 
properties of extended Lya nebulae ("Lya blobs") at 
z — 2.3. After searching over 4.82 deg^, we discover 
four Lya blobs with Llyq = 1-6 - 5.3 xlO^^ ergss"^, 
isophotal areas of 28 - 570", and broad Lya line profiles 
{Av = 900 - 1250 km s^^). We confirm the redshifts of 
all four blobs spectroscopically. In contrast with the ex- 
tended Lya halos associated with high-0 radio galaxies, 
none of our four blobs are radio-loud. The X-ray lumi- 
nosities and optical spectra of these blobs are diverse. 
Two blobs (3 and 4) are X-ray-detected with Lx(2— 7 
keV) = 2-4 xlO'''* ergss"-"^ and have broad optical emis- 
sion lines (C IV and C III]) characteristic of AGN, im- 
plying that 50% of our sample blobs are associated with 
strong AGN. The other 50% of blobs (1 and 2) are not 
X-ray or optically-detected as AGN down to similar lim- 
its, suggesting that AGN like those in blobs 3 and 4 are 
not necessary to power them. The number density of the 
four blobs is 3 x 10~^ hjQ Mpc~'^, comparable to that 
of galaxy clusters at similar redshifts and 3 x lower than 
that found in the SSA22 proto-cluster at z = 3.1, even 
after accounting for the over-density of that region. The 
two X-ray undetected blobs (1 and 2) are separated by 
only 70" (550 physical kpc) and have almost identical 
redshifts (corresponding to < 360 physical kpc along the 
line-of-sight), suggesting that they are part of the same 
system. Given the rarity of the blobs and our discovery 
of a close pair, we speculate that blobs occupy the high- 
est density regions and thus may be precursors of today's 
rich cluster galaxies. 
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TABLE 1 
Properties of Lva Blobs 



o 



ID 


R.A. 


Dec. 


z 


L(Lya) 


Size 


/recv 




L(2-7keV)'' 


NDWFS Name 


NDWFS magnitudes'' 




(J2000) 


(J2000) 




(10*^crgss-l) 


(□") 




(A) 


(lO'^'crgss-i) 




Bw 


R 


I K 


Blob 1 


14 30 59.0 


+35 33 24.7 


2.3186 


2.70 ± 0.10 


43 ± 4.8 


0.93 


422 ± 20.8 


< 0.63^= 


NDWFS J143059.0+353324 


23.63 


22.91 


22.69 


Blob 2 


14 30 57.8 


+35 34 31.6 


2.3178 


1.61 ± 0.08 


29 ± 6.9 


0.57 


383 ±22.2 


< 0.18': 


NDWFS J143057.8+353431 


24.05 


22.92 


22.53 


Blob 3 


14 37 25.1 


+35 10 48.8 


2.3321 


5.25 ± 0.19 


57 ± 5.2 


0.97 


171 ± 7.3 


4.11 ± 1.86 


NDWFS J143725.0+351048 


22.11 


20.67 


19.93 16.74 


Blob 4 


14 25 03.5 


+34 58 55.1 


2.3211 


2.27 ± 0.12 


28 ± 4.6 


0.67 


174 ± 10.4 


1.63 ± 1.40 


NDWFS J142503.4+345854<i 


23.06 


23.70 


22.13 




















NDWFS J142503.3+345855'^ 




23.66 


23.24 


X-ray luminosities are derived from the observed flux in the 0.5-2 keV band assuming a power-law spectrum with T = 


1.7 and a Galactic absorption ATh 


= 1 X 10 


2°cm-i. 





^ Vega magnitudes from the NDWFS DR3 catalog. 

The upper Hmits are derived from the combined Chandra ACIS-I images (Obs-ID 3608,6993). 
^ It is not clear which of two galaxies is associated with the blob 4 and the X-ray source. 
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